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We use the three-dimensional Josephson junction array system as a model for studying the tem- 
perature dependence of the c-axis resistivity of high temperature superconductors, in the presence of 
an external magnetic field H applied in the c-direction. We show that the temperature at which the 
dissipation becomes different from zero corresponds to a percolation transition of the vortex lattice. 
In addition, the qualitative features of the resistivity vs. temperature curves close to the transition 
are obtained starting from the geometrical configurations of the vortices. The results apply to the 
\ cases H =^ and H — 0. 
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1 , Strong thermal fluctuations and anisotropy make the 
~—> ' physics of the vortex lattice in high-Tc materials be much 
, more rich and complicated than predictions of mean field 
C*~) ■ theories Q. This shows up, in particular, in the com- 
plicated structure of the field-temperature (H-T) phase 
• diagram of the high-Tc's. It seems to be clear that 
, there is a line in the H-T phase diagram that sepa- 

■ rates a low-temperature phase (known as the vortex glass 
phase Q) where vortex lines are frozen in space, and a 

■ high-temperature phase in which the vortex lines move 
, through the material due to thermal activation. The 

■ passage to the normal state when increasing tempera- 
ture is likely to be a crossover, instead of a well-defined 
transition. The curve that separates the low- and high- 

, temperature phases is named irreversibility line (IL) . The 
^ ^ 1 V-I characteristics when an external current is applied 
O 1 perpendicularly to the magnetic field is different above 
and below the IL. Below IL the V-I curves are well fit- 
ted by V ~ exp [(— J C /J) M ] with p, and I c being two pa- 
rameters J^,||, and in particular, the resistivity p of the 
system -which is defined as p = limj_>o (V/ J)- i s strictly 
zero. Above IL the behavior is ohmic, i.e., V ~ J 

When the current is applied parallel to the field, the 
mean force exerted on the vortices is zero. However there 
are local forces -due to misalignment of the local mag- 
netic field- that may give rise to dissipation. The most 
important mechanism for dissipation in this configuration 
at intermediate temperatures is the thermal activation of 
vortex loops, which gives a voltage V c ~ exp [—I c /I] [Q, 
implying zero resistivity. In this work we show that when 
temperature is increased there is a phase transition at a 
temperature T p that reflects a thermodynamic property 
of the vortex system and is signed by the occurrence of a 
non-zero resistivity. In fact, the I — V characteristic for 
YBaCuO, when current and magnetic field are parallel 
to the c-axis show the following behavior H|| : for small 
currents and high temperature the response is ohmic, the 
range of currents that gives a linear response is reduced 
as the temperature decreases, and at a well-defined tern- 
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perature T p the linear behavior disappears. Moreover, 
the I — V curves can be scaled on two universal curves, 
corresponding to T > T p and T < T p respectively 0. 
This behavior -that is similar to what occurs when the 
current is applied in the a6-plane- supports the idea of a 
thermodynamic transition that we identify with a perco- 
lation transition of the vortex system. Experimentally, it 
is observed H that the dissipation in the c-axis appears 
at different temperatures than in the a6-plane. This im- 
plies that the 'irreversibility line' for a current parallel 
to the field is different than the corresponding to the ab- 
plane. 

Here we explore the following idea. At zero tempera- 
ture, the vortices are straight lines, and the net force on 
each of them when a small current in the c-axis is present 
is zero. At low temperature, vortex lines start to wander 
and vortex loops are created due to thermal activation. 
However, if temperature is not too high, vortex loops and 
vortex lines arc still isolated from each other and dissipa- 
tion in the linear regime is zero -except for surface effects 
(see below). When increasing temperature, vortex lines 
and thermally generated vortex loops start to touch each 
other and for temperatures greater than a critical value 
T p , there will be a vortex path crossing the sample along 
the a6-plane. The net force exerted by the current on 
this path is different from zero, and a finite dissipation 
will be observed. In this way we qualitatively see that 
the existence of paths perpendicular to the current in 
the sample -i.e., the transversal percolation of the vortex 
lattice- is crucial for the dissipation in the c-axis S . 

The model used to test this idea is the three dimen- 
sional Josephson junction array on a discrete lattice, that 
has been described in detail elsewhere PJnj ■ The dynam- 
ics of the model is contained in the evolution of the phases 
<p l (t), which are defined on the nodes of a cubic lattice 
and represent the phase of the order parameter. Between 
nearest neighbors nodes there are Josephson junctions 
characterized by a critical current Iq and a normal resis- 
tance Rq. The equations describing the model arc 

f = iosin (<p< - / - A«) + if-^^ 1 + rfHt) 

(1) 

■>'■"■ ( 2 ) 

{*'} 

Eq. |l| gives the current j u between nearest neighbors 
nodes i and i' with phases (p l and ip l . Here A n is the 
vector potential of the external magnetic field, and rf l (t) 
is an uncorrelated gaussian noise which incorporates the 
effect of temperature. Eq. ^| assures the current con- 
servation on each node, and j l ext is the external current 
applied at node i. 

The model allows for the existence of vortices, which 
consist in singularities of the phases ip(t) around a given 
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closed path. Self-inductance and disorder effects are not 
considered and the system is taken isotropic for simplicity 
-i.e., Iq and Rq are taken constant throughout the lattice. 

We numerically integrate Eqs. |^J| in time. Voltages 
at different points of the sample are calculated as the 
temporal mean value of the time derivative of the phases 
if. The resistivity of the sample in a given direction is 
calculated by injecting a small external current (typically 
around ~ 1/20 of the critical current of the junctions) by 
one of the faces of the sample and withdrawing it from 
the opposite face. The small value of the external current 
is chosen in order to be in the linear regime, in which the 
voltage drop is proportional to the applied current. 

The boundary conditions (BC) are taken open in the 
a&-plane. However, if open BC in the c-axis are used, 
there will be a finite force on an isolated vortex at finite 
temperature if the top and bottom ends of the vortex 
are not aligned. The dissipation -that is non-zero even 
in the linear regime- caused by this net force turns out 
to be independent of the thickness of the sample [fl0|| , 
and in this sense, it is only a surface effect. In order 
to eliminate this spurious surface effect it is crucial to 
use BC for the c-direction that assure that each vortex 
line leaving the sample at a given point of the bottom 
plane re-enters at the same point of the top plane. Strict 
periodic BC on the phases tp have this property, however 
we would obtain that the voltage difference between top 
and bottom planes is identically zero. We use, instead, 
open BC for the mean value of the phases in the top 
(<^r)and bottom (</?~b) planes, and periodic BC for all the 
phase differences tp l T — <pr T and ip B — Lp 3 B . This guarantees 
the periodicity of the vortex configurations and permits 
the calculation of the c-axis resistivity. 

We have to define a criterium for percolation: in our 
model there is a typical length which is the lattice pa- 
rameter a. Distances smaller than a cannot be resolved. 
Flux conservation implies that every flux line going into 
a unit cell of our lattice also goes out of the cell. When 
two vortices go into the same elemental cell we cannot 
tell which one of the two outgoing vortices correspond 
to each one of the ingoing vortices. We interpret this 
situation as the meeting of two vortex lines. In a real 
material this corresponds to two vortex lines being at a 
distance lower than the core size of the vortex. At high 
enough temperatures the vortex structure may percolate 
perpendicularly to the applied field: starting from one 
side of the sample we can follow a vortex line and arrive 
to the opposite side of the sample. Due to the finite size 
of the systems used, and to the dynamical evolution, per- 
colation is not expected to occur at every time, but only 
at a given fraction of the total time, which depends on 
temperature. We evaluate the probability that there ex- 
ists a vortex line crossing the system from one side to the 
opposite as a function of temperature. Because a sharp 
percolation transition can only be seen in the thermody- 
namic limit, we do scaling with the size of the system. 
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In Fig.[l](a) we show the resistivity of a cubic (L a b x 
L a bxL c , L a i, — L c = L) sample for an external field of 0.2 
(in units of quantum fluxes per plaquette) as a function 
of temperature (which is measured in units of the Joseph- 
son energy of the junctions) for three different sizes of the 
system: L = 8, 16 and 24. For comparison, the resistivity 
when the current is applied perpendicularly to the field is 
also shown for the case L = 8. It is clearly seen that the 
onset temperature for the dissipation in the c-axis Tp is 
higher than the corresponding to the a&-plane. Fig.Wb) 
shows the probability that the vortex lines have perco- 
lated through the sample along the afr-plane. We see a 
percolation transition around T p that becomes narrower 
the greater the size of the system. This indicates that 
there exists a sharp percolation transition in the ther- 
modynamic limit. As an additional check, in Fig. [j](b) 
(inset) the data of Fig. |l](b) are plotted vs a rescaled 

variable x: x = L a ah 1/2 - (1 - exp (-A/T)) ic , where 

a = 0.7, and A = 3.75 are numerically found parame- 
ters. This scaling comes up by using a simple model for 
the percolation p~T[ j . It strongly suggest that a (percola- 
tive) thermodynamical phase transition is occurring in 
the system. 

By comparing Figs. [l](a) and (b) it can be seen that 
the temperature T p where c-axis resistivity starts to be 
different from zero is the same temperature at which the 
percolation probability becomes finite. This indicates -as 
qualitatively discussed above- that the percolation tran- 
sition is a necessary condition for the existence of dissi- 
pation in the c-direction. 

In addition, we would like to have a more quantitative 
estimation of the resistivity, based on the geometrical 
configurations of the vortex system. This can be accom- 
plished in the following way: Let us consider a sample 
of size L c (L a b) in the c(a6)-direction The resistivity p 
of the sample in the c-direction is proportional to the 
number of paths n that cross the sample in the afo-plane 
per unit of area, times the velocity v this paths acquire 
under the external force, divided by the external current 
density j: p ~ nv/j. The velocity v is given -using a vis- 
cous fluid argument- by the external force F divided by 
a total viscosity 77, which is equal to a specific viscosity 
coefficient, ijo times the total length of the vortex path, 
that we will call I, i.e., v — F/rjol. The force F is given in 
term of the external current and the size of the system: 
F ~ jLat,. We obtain p ~ nL a {>/riol. The coefficient 770 
depends on temperature, however, on small ranges near 
the percolation threshold we will take it as a constant. 
The determination of n and I is a difficult task, because 
the percolation paths across the sample are not uniquely 
defined due to the crossing of vortex lines (see Fig.|(a)). 
We will use the following estimation: we assume that 
nxl x L a b x L c is the volume S of the percolation cluster 
in a sample of volume L Q b x L a b x L c . The value of S can 
be easily evaluated from the numerical simulation. We 
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obtain p ~ S/r]ol 2 L c . It remains to estimate the value 
of I. This length I depends both on temperature and 
the size of the system. As we said, a direct numerical 
determination of I is difficult due to indeterminacies at 
the crossing points of the vortex lines. We will use the 
most crude estimation (see Fig.||): when the magnetic 
field H is close to zero -i.e. H < H cross , where H cross is 
a crossover field which is defined below- we take I ~ L ao . 
However, for H > H cross , percolation proceeds via the 
external field generated vortices and the length of a per- 
colation path is much larger, and can be estimated to be 
I ~ LcLab/H- 1 / 2 . In this way we obtain the following 
scaling for the resistivity near the percolation threshold 

p~S/L 2 ab L c for H<H cross , (3) 

p~S/L 2 ab Ll for H>H cross . (4) 

This scaling is expected to be valid only close to the 
percolation threshold. 

The crossover field H cross is estimated as H cross ~ 
and corresponds to the zero-temperature lattice 
parameter of the vortex structure being equal to the 
thickness of the sample. (The value of this field is about 
20 G for a 1 fim thick sample). For the value H = 0.2 
used in Fig. [l] we are in the case H H cross for all the 
values of L z considered. 

In order to check the previous estimations, in Figj^(a) 
we compare the values of pL 2 ab h\ and S vs temperature 
when L ab is varied between 16 and 30, for H — 0.2. In 
Fig. ||(b) pL 2 ah L\ and S vs temperature are compared 
when L z is varied between 12 and 24, for the same field 
H = 0.2. The only free parameter of the fitting is a 
global factor, which is the same in Figs. |(a) and |(b). 
The agreement between the numerically calculated val- 
ues and the estimated ones close to the threshold is fairly 
good if we take into account all the approximations made 
in order to obtain Eqs. || and |J. A more precise estima- 
tion of the resistivity using only the geometrical configu- 
rations of vortex lines seems to be difficult because of the 
following facts: The percolation paths across the sample 
are not uniquely defined (see Fig. ||(a)), and the real 
movement of vortex lines under the external force will 
depend on the cutting energy. The viscosity i]q is not 
a constant, but a function of temperature. In addition, 
the supposition of a phenomenological viscous motion of 
vortex lines may not be accurate at low temperatures, 
when vortices creep. 

The existence of two resistive transitions (in the c-axis 
and the aft-plane) has been experimentally observed in 
YBaCuO [pi. The values of the two characteristic tem- 
peratures depend on the pinning, vortex elasticity and 
magnetic field. In YBaCuO, as the thickness of the sam- 
ple increases the two temperatures become closer to each 
other. In our simulations we find that the temperature 
at which the percolation transition occurs decreases as 
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~ 1/ ln(L c ) jy]], as it can be deduced from the scaling in 
Fig. 0(b) (inset). 

The thermal excitations in the form of vortex lines 
crossing the sample along the aft-plane destroy the phase 
coherence along the c-axis. For T > T p the coherence 
length £ c is of the order of the mean distance between per- 
colation paths, i.e., £ c ~ Lc/n 1 / 2 . We conclude that the 
mechanism that leads to the 2D-3D transition in high-T c 
materials with moderate anisotropy is the percolation of 
vortex line perpendicular to the external field. 

In summary, for a model high temperature supercon- 
ductor we have shown by using qualitative arguments and 
numerical simulations, that the onset of the resistivity in 
the c-direction is related to a percolation transition of 
vortex lines in the aft-plane. The results hold for H ^ 
and H = 0. A qualitative estimation of the resistivity 
near the threshold, and its finite size scaling has been 
given. For the sizes of the isotropic systems used, the 
percolation transition occurs at higher temperature than 
the resistive transition in the aft-plane, and corresponds 
to a new thermodynamic transition that should be char- 
acterized by new critical exponents different from those 
obtained for the vortex glass transition when current is 
applied parallel to the aft-plane. We expect these results 
to be valid also for anisotropic systems, at least in the 
case of moderate anisotropy, as in YBaCuO. 
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FIG. 1. (a) Resistivity along the c-axis (solid circles) and 
along the ab- plane (open circles), and (b) probability of per- 
colation across the sample, for a cubic lattice of size L vs. 
temperature (in units of the Josephson energy of an individ- 
ual junction). Inset: percolation probability vs scaled temper- 
ature x (see the text for definition) Different symbols corre- 
spond to different sizes of the sample. 

FIG. 2. Schematic two-dimensional representation of per- 
colation paths (broken lines) of vortex lines (solid lines) across 
the sample for the cases H — (a) and H > (b). Note in (a) 
the two different percolation paths with the same end points, 
and in (b) the percolation through the externally generated 
vortex lines. 

FIG. 3. Numerically calculated resistivity times L ab L'l 
(solid symbols) and volume of the percolation cluster S (open 
symbols) vs temperature for different sizes of the sample 
(Lab x Lab x L c ), and H = 0.2. 
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